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Abstract. In this paper we present a multidimensional version of the van der 
Corput lemma where the decay of the oscillatory integral is gained with respect to 
all space variables, connecting the standard one-dimensional van der Corput lemma 
with the stationary phase method. 

1. Introduction 

This paper is devoted to the estimates for the oscillatory integrals of the type 

/(A) = f e iX * {x) a(x) dx , 



where the support of a e C™(M. N ) is sufficiently small. The estimate for /(A) as 
A — > oo is well-known in one dimension N — 1. If $ is real-valued and > 1 

on the support of a, then the estimate |/(A)| < c^A -1 ^ holds when k > 2, or when 
k = 1 and is monotonic. In this case the bound Ck is also independent of $ and 
A (see e.g. Sogge [TDJ or Stein [H]), and the decay rate is sharp. This result plays a 
crucial role in various areas of analysis. For example, it is closely related to sublevel 
set estimates of the form 

meas {s E supp a : |$(s)| < t} < c k t 1/k , 

where $ is a function as above, with numerous applications in partial differential 
equations, microlocal analysis, harmonic analysis, etc. 

A multidimensional version of these results would be of great value, but presents 
many difficulties. It is known that for dimensions N > 1, if, for example, <9 a $ > 1 
on supp a, then |/(A)| < CalA!" 1 /'"'. The decay rate here is sharp, but the constant c a 
may depend on $ and the estimate does not scale well. Again, such estimate is closely 
related to the multilinear sublevel set problem (see e.g. Phong, Stein and Sturm [7]). 
Parameter dependent sublevel set estimates were recently used by Kamotski and 
Ruzhansky [5] in the analysis of elliptic and hyperbolic systems with multiplicities, 
to yield Sobolev space estimates for relevant classes of oscillatory integrals and for 
the solutions. 

Recently, Christ, Carbery and Wright [2 J and Carbery and Wright [3], proposed 
versions of van der Corput lemma for functions of several variables, in formulations 
where the constant in the estimate is independent of the phase function. This aspect 
is of significant importance for applications allowing to investigate various pertur- 
bation and other properties of appearing integrals. However, the decay rate of the 
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corresponding oscillatory integral there is again essentially one-dimensional because 
the non-degeneracy of only one (higher-order) derivative is assumed. 

At the same time, the decay rate exhibited in many problems of interest is better 
than one- dimensional. If one compares this with the case of non-degenerate stationary 
points of $, the stationary phase method will readily yield the decay rate |/(A)| < 
(j\-N/2 However, if a stationary point degenerates, the situation becomes much more 
delicate (see e.g. Hormander [U Chapter 7]), and no good estimates are available in 
general. 

The aim of the present paper is to bridge the gap between van der Corput lemma 
and estimates provided by the stationary phase method. On one hand, the stan- 
dard van der Corput lemma works well for degeneracies of high orders but produces 
only one-dimensional decay rate. On the other hand, the stationary phase method 
produces the multidimensional decay rate, but does not work well for degenerate 
stationary points. 

The result of this paper will yield a multidimensional decay rate for degenerate 
stationary points. We will identify a class of functions, for which this can be achieved. 
These functions will have certain convexity type properties. It is clear that certain 
convexity conditions are necessary to ensure the multidimensional decay rate. In fact, 
conditions of the one-dimensional van der Corput lemma guarantee that the function 
(or some derivative of the function) is convex in one dimension. Thus, it is natural 
that an analogue of convexity also appears in several dimensions to ensure that we 
gain one-dimensional decays in all directions. It is then a question of putting all 
these rates together to yield the full multidimensional decay, which will turn out to 
be iV-times better than the standard van der Corput estimate. 

In what follows we will also allow phase function $ to be complex valued and to 
depend on an arbitrary set of parameters. These two situations often happen in ap- 
plications to partial differential equations, in particular in the analysis of solutions 
represented as oscillatory integrals, leading to the dispersive and to the subsequent 
Strichartz estimates. Thus, the complex phase corresponds to the fact that charac- 
teristics of the analysed evolution equations may be complex (see e.g. Treves [13]). 
At the same time, the dependence of the phase and of the amplitude on parameters is 
also essential, and is related to uniform estimates of [2J. Also, in applications to the 
Strichartz estimates for hyperbolic equations of high orders considered by Ruzhansky 
and Smith [9], a parameter is essential to encode the information on low order per- 
turbations of the equation, in order to establish the dispersive estimates for solutions 
uniformly over such perturbations. At the same time, in hyperbolic equations with 
time dependent coefficients (e.g. considered by Matsuyama and Ruzhansky [B]), the 
parameter encodes the information on the perturbations of the limiting behaviour of 
coefficients, again allowing to obtain dispersive estimates uniformly over such per- 
turbations. We will leave out these and other applications outside the scope of this 
short paper. 

We will use the standard multi- index notation a = (ai, . . . , ajy), denote its length 
by | a | = ai + • ■ ■ + a/v and partial derivatives by d a = d"l ■ ■ ■ d%£- We will also use 
the standard convention to denote all constants by letter C although they may have 
different values on different occasions. 



MULTIDIMENSIONAL DECAY IN VAN DER CORPUT LEMMA 



3 



2. MULTIDIMENTIONAL VAN DER CORPUT LEMMA 



The following theorem is the main result that establishes the multidimensional 
decay rate for a class of oscillatory integrals. 

Theorem 2.1. Consider the oscillatory integral 



I{\v) 



where N > 1, and v is a parameter. Let 7 > 2 be an integer. Assume that 

(Al) there exists a sufficiently small 5 > such that \ G (-85/2(0)), where 

-85/2(0) is the ball with radius 5/2 around 0; 
(A2) $(x, v) is a complex valued function such that Im$(x, v) > for all x G 

supp x an d all parameters v, 
(A3) for some fixed z G supp x, the function 

F(p, to, v) := <&{z + puj, u), \u\ — 1, 

satisfies the following conditions. Assume that for each p = (uo,v), func- 
tion F(-,p) is of class C 7+1 on supp x, an d let us write its j th order Taylor 
expansion in p at as 

7 

j=0 

where -R 7+ i is the remainder term. Assume that we have 
(Fl) a (p) = a\(p) = for all p; 

(F2) there exists a constant C > such that Yl]=2\ a j(h L )\ — ^ f or °^ /•*> 
(F3) /or eac/i p, \d p F(p, p)\ is increasing in p for < p < 5; 
(F4) /or each k < 7 + 1, dpF(p, p) is bounded uniformly in < p < 5 and p; 
(A4) /or eac/i multi-index a of length \a\ < [^] + 1, t/iere exists a constant C a > 
suc/i t/iat |<9"a(x, z/)| < C Q /or a// x G supp % and all parameters v . 
Then there exists a constant C = Cjv, 7 > such that 

— — 

(2.1) I ( A, u)\ < C(l + A) t /or a// A G [0, 00) and all parameters v. 

Theorem 12.11 obviously includes the case where a and $ depend on different sets of 
parameters. In this case we may let v run over the whole space of parameters. 



We also note that assumption (A3), or rather (F3) , can be view as an analogue of a 
convexity assumption. Indeed, if F is real valued, then (F3) implies that the second 
order derivative d^F(p,p) does not change sign for < p < 6, because d p F(0,p) = 
by (Fl). In turn, condition (Fl) is not restrictive, since ao(p) can be taken out of 



the integral, and non-zero a\(p) would actually give a faster decay rate. 

Proof. It is clear that (12. II) holds for < A < 1 since |/(A, v)\ is bounded for such 
A, in view of assumptions (Al), (A2) and (A4). So, we may consider the case where 
A > 1. Let z G M. N be as in (A3), and set x = z + poo, where uj G S^ -1 , p > 0. For 
A = 1, we use §° = {—1, 1}. Then we can write 



^ t+pu ^a(z + pu, u) X (z + poo)p N - 1 dp du . 
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It suffices to prove ( 12. 1ft for the inner integral. 

Choose a function 9 G Cq°([0, oo)), < 9(s) < 1 for all s, such that 0(s) is 
identically 1 for < s < | and is identically zero for s > 1. Then with our notation 
F(p, u, v) = $(z + poo, u), we split the inner integral into the sum of the two integrals 



I 2 (\,u,uj,z) = 







D i\F(p,u),i/) 



a(z + puj,v)x(z + pu)(l - 9)(\~< p)p dp. 



Let us first estimate I± = ii(A, v, u>, z). Since 9(\f p) = for \i p > 1, changing 

x 

variable r = A^p, we have 



\h\<C J 9(\^p)p N - L dp = C I t"-^-* 
which yields the following estimate for J^: 



(r)A t dr, 



(2.2) 



In order to estimate I 2 = ^(A, v, uj, z), let us first establish a useful estimate for 
functions F satisfying condition (F3) We claim that under condition (A3) , or rather 
under (F1)-(F4), there exist constants C,C m > such that we have estimates 

\d p F{p,n)\>Cp^ 



(2.3) 
(2.4) 



and \d™F{p,p)\ < Cnf^dpFip, n)\ 



for all < p < 5, all parameters p, and all m < 7 + 1. First, note that for < p < 1 
and m — 7 + 1, estimate (12.41) follows from (I2.3P and assumption (F4) So we may 
only consider m < 7. 

Now, assumption (F2) implies that 

(2.5) 7r(p,p) := > ■ 

J=2 

Thus, in order to prove (12.31) . it suffices to show that 

(2.6) \d p F(p, p)\ > Cir(p,p) for all < p < 5 and all parameters p. 



For 1 < m < 7, we have, using (A3) 
(2.7) 



7— m 



9™F(p, P) = J2 fcj Q fc+m(/^)p fc + #m, 7 -m(p, M) 



fc=0 

■V :YY+1 cV„ ■■\(p-a) 7 - m 



where R mn - m (p, p) = f£ d] +1 F(s, p)^j^^jr ds is the remainder term of the (7— m) 
Taylor expansion of d™F(p,p). By (F4) and (\2.5L we get that 



th 



(2.1 



|iC, 7 - m (p,/i)| < C 7 , m p^ m+1 < C 7im vr(p,p)p 2 - m for < p < 5. 
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Hence, for < p < 5, we have 

7-1 

\d p F(p,p)\ = ly^ik + l)a k+l (p)p k + i2 li7 _!(p,/x) 

fc=0 

i 7 i 7 



i=2 



i=2 



It follows now from assumptions (Fl) and (F3) that 
\d p F{p,fi) \ = max \d p F{a,p)\ 

0<cr<p 



> 



max VS jaAp)(ji 1 — max C 7 7r(er, p)a 



0<CT<p 



i=2 

7 



\)<~o<~ 1 1 * 

J=2 



since 7r(cr, yu)cr = X!J=2.?'l a i(/ i )l cr ' ? achieves its maximum on < a < p at cr = p. 
Noting that 



tax > 

i=2 



max 

o< 



and 

i=2 

are both norms on C 7_1 and, hence, are equivalent, we immediately get 



\d p F{p,fi)\ >C^2j\aj{fi)\p3 1 - C 1 ir(p,p)p 

3=2 

>(C-C 7 5)7r(p,//)>C7r(p,//), 

for some constants C > 0, if 5 is sufficiently small. This completes the proof of (12.61) . 

To prove f)2.4p . we will use the representation (I2.7p . Since 1 < m < 7, it follows 
from the definition of 7r(p, p) that 



fc=0 



which, together with (12. 8ft and ( 12. 6ft . yields 

< C^y-^ia^Op,^)! for < p < 5. 

This completes the proof of the claimed estimates (12. 3ft and ( 12.41) . 
Let us now come back to the estimate for I 2 . Define the operator 

L:= {i\d p F{p,u,v))- x ^ 

which clearly satisfies the useful identity L(e* AF ( p ' w '^) = e 1 ^^'^ . Denoting the 
adjoint of L by L*, we have, for each IgNU {0}, 



i\F(p,ui,u) 1 t *\l 



(L*) l [a(z + pu, v) X {z + pu)(l - 6)(\^p)p N - 1 ] dp . 
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Now, 



(L 



,s p ,p,r,l ' 



d Sl F- 



d p p F 



(dpF) l +P 



where the sum is over all integers s\ 
From (E3D and (E3D it follows that 

d s /F 



SpiP, r > such that Si + - ■ ■ + s p + r —p = I. 



d?F. 



(d p F) 



i+p 



(j>,u,u) <Cff 



9 — Sl S p —l~f+l 



Cp 



r—l-y 



Also, it is easy to see that for r < [— ] + 1, we have 



(2.9) 



dp~' 



ia(z + pu J ,v)x(z + pu J )(l-6)(\^p)p"- L } < Cnp^-^p) 



where x(A,p) is a smooth function in p which is zero for \i p < |. Let us now take 



1 = 1 



1, so that iV — £7 < 0. Then we can estimate 

POO 

\h\ <cva ' / Csi,...,s p ,p, 

Jo 

J^-^dp 



;i P' 



^p N ~ l ~r X{\p)dp 



<C N X 



A T 



CmX 



p 



N-l-y , 



U 7 



c 



AT, 7 



JV 



Combining this estimate with estimate (12.21) for Ii, we obtain the desired esti- 
mate ( 12. ip . This completes the proof of the theorem. □ 



We note that in the proof we showed that if function F satisfies conditions (Fl) 



(F4) , it also satisfies estimates ( 12.31) and (12.41) . A version of this part of the argument 



was discussed by Sugimoto [12] for real valued analytic functions without dependence 
on p, where the analysis was based on the Cauchy's integral formula for analytic 
functions (see also Randol [8] and Beals pQ). The proof that we give for ( 12.31) and 
(12. 4p extends it to the generality required for Theorem 12.11 

In fact, let us also briefly indicate a smooth version of these estimates. Suppose 
that a function F(-,p) is smooth in the first variable, and that it satisfies conditions 
(Fl) - (F3) , as well as condition (F4)| for all m £ N. Then we claim that for sufficiently 



small 5 > 0, estimates ( 12.31) and ( 12. 4p are satisfied also for all m £ N. 

Indeed, we already proved estimate ( 12. 3D and we also proved ( 12.41) for m < 7. It 
remains to consider the case m > 7. Since 7 + 1 — m < 0, from (F4) it trivially 
follows that for < p < 5 we have a stronger estimate 



\d™F(p,p)\<C m <C m ^ +1 - m 
where the last inequality is a consequence of (12.31) . 



<C mj5 p'~ m \dpF(p,p)\, 
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